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Abstract 

Understanding the ramifications of reduced crystalline symmetry on magnetic behavior is 
a critical step in improving our understanding of nanoscale and interfacial magnetism. How¬ 
ever, investigations of such effects are often controversial largely due to the challenges inherent 
in directly correlating nanoscale stoichiometry and structure to magnetic behavior. Here, we 
describe how to use Transmission Electron Microscope (TEM) to obtain Electron Magnetic 
Circular Dichroism (EMCD) signals as a function of scattering angle to locally probe the mag¬ 
netic behavior of thin oxide layers grown on an Fe (1 1 0) surface. Experiments and simulations 
both reveal a strong dependence of the magnetic orbital to spin ratio on its scattering vector 
in reciprocal space. We exploit this variation to extract the magnetic properties of the oxide 
cladding layer, showing that it locally may exhibit an enhanced orbital to spin moment ratio. 

This hnding is supported here by both spatially and angularly resolved EMCD measurements, 
opening up the way for compelling investigations into how magnetic properties are affected by 
nanoscale features. 

Systems of restricted size and dimensionality represent a frontier for research on magnetic ma¬ 
terials. By reducing the spatial dimensions of crystalline magnets, interfacial magnetic properties 
become more prominent in the overall magnetic behavior of the system^. Controlling interfaces 
through the fabrication of magnetic heterostructures enables researchers to produce materials that 
exhibit entirely new properties from the bulk constituents alone. For these reasons, considerable 
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research efforts are currently dedicated to improving the understanding of nanoscale magnetic be¬ 
havior; however, many of the techniques capable of quantifying magnetic moments - such as X-ray 
Magnetic Circular Dichroism (XMCD) - lack the spatial resolution necessary to directly correlate 
the magnetic behavior of the material to the nanoscale features from which they arise. 

An example of a field where such understanding is necessary is research on magnetic transition 
metal oxides, particularly the iron oxides. Iron oxides form in a variety of phases exhibiting a 
wide range of magnetic behavior. Of particular interest is the phase Fe 304 , commonly known as 
magnetite. Magnetite has potential to play a role in the development of nanoscale magnetic appli¬ 
cations due to its anticipated half-metallic behavior® leading to a near 100% spin polarization®!^ 
high chemical stability in ambient conditions, and ability to stabilize very thin films against the 
onset of superparamagnetism®. It crystallizes in an inverse spinel cubic structure with Fe assuming 
both +2 and -|-3 oxidation states in a ferrimagnetic arrangement. Fe^"*" is evenly distributed among 
the tetragonal and octahedral sites and is aligned antiparallel while Fe^"*" is found exclusively on 
the octahedral site, resulting in a net magnetic moment of nearly 4 /ib. It is generally considered 
that the orbital component of this net magnetic moment in bulk magnetite is nearly completely 
compensated by this symmetry, resulting in a very small net contribution to the total magnetiza¬ 
tion®®. However, it has been suggested that local structural and chemical variations may break this 
compensation, increasing its detectable magnitude®^. Since these effects are necessarily nanoscale 
in origin while the measurement techniques employed probe macroscale materials, it has proven 
exceedingly difficult to explore these effects in greater detail. 

One technique capable of contributing to this effort is Electron Magnetic Circular Dichroism 
(EMCD) in the Tranmission Electron Microscope (TEM). First proposed in 2003^ and experi¬ 
mentally demonstrated in 2006®, EMCD is an exceptionally useful technique for quantifying the 
magnetic properties of materials on the nanometer scale®® and has been employed in the analysis 
of magnetic domain walls at the nanoscale®. Magnetospirillum magnetotacticum®, LaSr-2 x 4 
manganese oxide nanowires®, Cr02 thin films®, Fe 304 nanopart ides®, and FeCo alloyd®. The 
technique uses the TEM specimen as a beam splitter in the electron microscope where, in the cor¬ 
rect scattering geometry, the fast electrons interact with the sample in a similar way to circularly 
polarized x-rays, enabling many of the measurements possible in XMCD to be executed in the 
TEM. Whereas an XMCD signal is obtained from two spectra acquired using right and left-handed 
circularly polarized beams, in EMCD, two Electron Energy Loss Spectra (EELS) are acquired at 
two conjugated scattering vectors. The EMCD signal itself is the difference between these two 
spectra. 

While the high spatial resolution of the EMCD technique clearly distinguishes it from XMCD, 
there are a number of additional differences that are just beginning to be explored and understood. 
First, since the EMCD signal depends on the electron channeling conditions of the sample in the 
TEM, it is possible to use this technique to obtain site-specific magnetic information ® — I Second, 
it has been demonstrated that the magnitude of the EMCD signal (measured as its signal to noise 
ratio) varies as a function of scattering angle (or q-vector)®®. This dependence - known as an 
“EMCD strength map” - can be experimentally recorded in the TEM either by angularly selecting 
the inelastically scattered electrons or by energy filtering a series of electron diffraction patterns. 
When sum rules are applied to the resultant signals, it becomes possible to quantitatively extract 
the orbital to spin magnetic moment ratio mi/mg as a function of q-space®®®. Significantly, 
while the magnitude of the EMCD signal changes as a function of crystalline symmetry of the 
magnetic material, the value of rriL/ms remains invariant for homogeneous systems®® 

In this communication, we exploit both of these dynamical diffraction properties to extract 
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for two magnetic materials having the same magnetic species but different crystal structures 
that overlap in the direction of the electron beam. Such situations readily arise for metallic TEM 
lamellae prepared in cross section that have been exposed to atmosphere between sample preparation 
and transfer to the microscope, forming a thin metal oxide layer on the exposed surfaces. If both the 
film as well as its oxide have a net magnetic moment, two independent, dissimilar EMCD strength 
maps will be generated. To a good approximation, the total EMCD signal detected at any given 
q-vector will thus be a linear combination of these EMCD signal maps. Since both EMCD strength 
maps vary differently in reciprocal space, when sum rules are applied to the entire EMCD signal 
map, the measured value for mi/ms will vary as a function of detector position, which can be 
experimentally determined with high precision. 

We exploit this effect to probe a TEM lamella of bcc iron having an exposed (1 1 0) surface, upon 
which a thin layer of cubic iron oxide has grown in the {111} orientation. Based on structural and 
spectroscopy investigations, we argue that this layer is best described as Fe 3_504 on both exposed 
surfaces, where 5 varies between 0 (yielding mixed valence magnetite, Fe 304 ) and 0.33 (yielding 
monovalent maghemite, 7 -Fe 203 )®. By selectively probing different regions in reciprocal space, we 
show that it is possible to extract magnetic information pertaining to both the underlying iron 
him as well as its thin oxide surface layer. Our experimental data combined with simulations on 
the system bring us to the conclusion that these thin Fe 3_504 layers may locally exhibit a large, 
uncompensated orbital magnetic moment. 


Results 

Analysis of iron surface oxidation 

A him of bcc Fe with a thickness of 50 nm was epitaxially grown on a single crystalline MgO 
substrate as detailed in the methods section. An initial quality assessment of the him using x-ray 
dihraction techniques verihed a close to single crystal (0 0 1) growth with a Full-Width at Half- 
Max (FWHM) of the rocking curve equal to 0.5°. A thin cross-sectional lamella was subsequently 
prepared for the TEM using the FIB in-situ lift-out method explained in the methods section. An 
overview and structural assessment of the TEM-prepared sample is provided in the supplementary 
section. 

In the thinnest area of the lamella, a region of interest was sought out for the EMCD measure¬ 
ments. A survey image of this region acquired with the microscope in Scanning TEM (STEM) mode 
using the High Angle Annular Dark Field (HAADF) detector is shown in figure The region used 
for this investigation is depicted in green. The elemental composition of the iron cross-section in 
this region was investigated with on-axis EELS measurements as described in the methods section. 
In the lower panel of figure real-space maps revealing the relative percentage of both Fe (rpe) and 
O (ro) are presented along with the calculated total absolute thickness (tabs) of the film for any 
given pixel position. The thickness of both the oxide cladding layers (tox) as well as the underlying 
metallic layer (tpe) was calculated under the assumption of the oxide being structurally similar to 
Fe 304 . See the methods section for details about the thickness calculation. 

The on-axis EELS maps reveal that the entire Fe film within the scanned region exhibits a 
detectable oxygen signal. In the middle of the Fe him, close to the substrate but slightly offset 
from the interface, oxygen comprises 10 - 20% of the total atomic concentration. This rises to 
nearly 60% in the upper region of the him, close to the vacuum. The absolute thickness of the 
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film decreases as a function of distance from the MgO substrate, which is consistent with a wedge 
shape, as expected from the sample preparation procedure. We therefore conclude that the oxide 
layer has encapsulated the original iron him and, towards the surface, has consumed nearly all of 
the metallic iron. 

Each pixel position in hgure [T] contains both an individual EELS spectrum as well as spatial 
information with a resolution of approximately 1.5 nm. We note that spreading of the beam with 
convergence angle a = 1.6 mrad during the propagation through a 50 nm sample thick is only 
approximately 50 nm ■ a = 50 nm ■ 0.0016 <0.1 nm. This property allows for the EELS data to be 
spatially segregated into regions corresponding to signihcant local features. Since the EMCD data 
were acquired from the same scanning region as the data in hgure [TJ these core-loss EELS data can 
be used to determine the amount of oxide that contributes to the corresponding EMCD signal. We 
thus identify two core regions of interest for the EMCD analysis. The hrst is a summation over all 
spectra in each spectral image. The second is a summation over the portion of the him with the 
highest oxide content, denoted the “surface” region. The area of summation is presented in hgure 
□ 

The structure of the oxide-encapsulated Fe layer was investigated with both Convergent Beam 
Electron Dihraction (CBED) and High Resolution TEM (HRTEM) techniques. The CBED pattern 
is shown in hgure This pattern was extracted from the middle of the boxed region shown in hgure 
[T] and thus reveals the structure of same area used for the EMCD measurements. The most intense 
rehections can be indexed as Fe [110], and this is done in hgure [^. The additional rehections 
appear to come in sets of three overlapping discs. The midpoint of these sets of rehections matches 
well with the indices expected from stoichiometric Fe 304 [111], and these labels are used to index 
the structure in hgure |^. The satellite rehections can be understood as a consequence of dynamical 
dihraction ehects resulting from the propagation of an electron wave through multiple crystalline 
lattices. This produces an interference lattice that manifests itself as superposition of an additional 
spatial frequency over all of the primary Bragg rehections. The frequency of this interference lattice 
Ag is 0.98 nm“^ corresponding to giio(Fe) — g422(FG304)) ^ind its lower orders are indexed in hgure 

Efc. 

To further rehne the structural analysis. High Resolution TEM (HRTEM) experiments were 
performed on a neighboring region of the him, as displayed in hgure The lower convergence angle 
ohered by this technique allows for sharper spots in the Fourier transform of HRTEM images than 
achievable in the CBED pattern. In hgure [^, the MgO substrate appears to be mono crystalline 
but the iron thin him is not. The discrete Fast Fourier Transform (EFT) of this image is shown 
in hgure |^, revealing that both high and low frequencies are present. The lowest order observed 
spatial frequency Ag is 0.98 nm“^, and this is visible in hgureas a low frequency “beat” oriented 
perpendicular to the substrate. This Ag vector can be attributed to an interference phenomenon 
known as Moire contrast. The analysis of this phenomena in this manuscript follows the approach 
described by Amidroi®. Critically, this approach demonstrates that the best explanation for the 
oxide structure is the mixed valence Fe 3_504 where S is close to 0. If 6 were close to 0.33 yielding 
a structure more similar to monovalent 7 -Fe 203 or maghemite, then one would expect additional 
rows of interference rehections due to the reduction of symmetry arising from the presence of 
cation vacancies needed to maintain charge neutrality®. Other iron oxides can be ruled out due to 
lacking the appropriate lattice plane spacings and symmetry operations to ht with the experimental 
observations. A more detailed analysis including simulations is provided in the SUPPLEMENTARY 
INFORMATION. We cannot exclude the possibility that minute amounts of other oxide structures 
or magnetite with impurities or vacancies may coexist with the Fe 3_504 in our sample. However, 


4 


these must be present in minute amounts so that they remain undetectable in the HRTEM images 
and diffraction patterns. As only large volumes of textured crystals can strongly modify the EMCD 
signal, we would thus expect the influence of any such phases to be negligible. 

Thus, based on the EELS, CBED, and HRTEM data, it can be concluded that the lamella is 
best described as a free-standing trilayer of Fe cladded between two thin Fe 3_504 layers where 6 
is close to 0. These layers appear to grow with a well-defined texture where Fe [1 1 0] || Fe 3_504 
[1 1 1]. R is worth noting that it has been previously shown that Fe 3_504 in this orientation is the 
most likely oxide to form on an Fe (1 1 0) surface under similar growth conditions^^^®. A schematic 
model of this composite structure is presented in figure |^. 

EMCD measurements 

The EMCD experiments were performed in a two-beam condition with g = Fe (0 0 2), as described 
in the methods section. A schematic of the diffraction pattern under these conditions is presented 
as an inset in figure |^,b. Two regions of q-space were sampled, and these regions are denoted as 
aperture pairings A and B throughout the text. Aperture pairing A was located on the Thales circle 
(figure]^ inset), while pairing B was located closer to the Fe (0 0 2) reflection (figure]^ inset). 
Each aperture pairing consisted of two individual aperture positions mirrored about the systematic 
row, thus collecting spectra corresponding to scattering geometries with opposite chirality. These 
positions are denoted Chiral Plus and Chiral Minus and are defined for each aperture pairing in the 
insets of hgure |^,b. 

Figure 1^ depicts the two chiral EELS spectra acquired using the scattering geometry of aperture 
pairing A, while figure shows the corresponding spectra collecting from aperture pairing B. All 
spectra in £gures|^ and b are a summation of all of the individual spectra in their respective spectral 
image and are plotted on a y-axis that has been normalized to the maximum value of the Chiral 
Plus spectrum in each pairing to facilitate comparison. In these figures, the pre-edge background for 
the Chiral Plus and Chiral Minus spectra has been subtracted and the signals have been shifted and 
aligned with respect to each other as described in the methods section. They otherwise represent 
the raw data. The difference between the two EELS spectra of opposite chirality for each aperture 
pairing is denominated as the dichroic signal (or “EMCD signal”) and is shown along with the 
spectra on the same normalized scale. The integral of both EMCD signals is presented in hgure 
A clear difference can be seen between the two aperture pairings, indicating that the ratio of 
the area under the Fe L 2 and Fe L 3 edges changes between aperture pairings. When sum rules are 
applied to the resulting EMCD signals, the value for mi/ms can be extracted as described in the 
methods section and the results for the two aperture pairings are shown in hgure |^. The values 
stabilize after approximately 745 eV and the remaining variations can most likely be attributed to 
the residual noise. The average values of mi/ms are 0.082 ± 0.008 for pairing A and 0.058 ± 0.004 
for pairing B. 

This variation of ms/ms in q-space strongly suggests that there are two sources contributing 
to the measured EMCD signal with dissimilar spin and orbital magnetic components that scatter 
diherently. The presence of two chemically and structurally distinct layers described in the previous 
section ohers a plausible explanation for this ehect. To more closely correlate this ehect to the 
diherent layers, we sum the spectra over the “surface” region denoted in hgure [T] The results of this 
summation for aperture pairings A and B are presented in hguresj^ and[^, respectively. A striking 
variation is visible. For aperture pairing A, the asymmetry on the Fe L 3 edge is approximately 10% 
while on the Fe L 2 edge this asymmetry is below the noise level. This will lead to an increased value 
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of mi/ms- This asymmetry nearly vanishes in aperture pairing B, which is additional evidence that 
the EMCD signal in aperture pairing B is dominated by the metallic iron signal. 

EMCD simulations 

To explore potential origins for the q-space dependence of mi/ms in this composite system, we 
simulated its inelastic electron scattering behavior. Figures and b show calculated maps of the 
magnetic signal originating from L 3 edge of iron, separately for the iron layer (fig. |^) and for the 
oxide layers (fig. g>). This signal corresponds to the difference spectrum between chiral minus and 
chiral plus positions of the L 3 edge, as described elsewhere^. It shows that the magnetic signal 
varies with scattering angles differently for iron than for the oxide layers. Note also the dissimilar 
order of magnitude of the signals resulting from the different thicknesses and densities of iron atoms 
between the two layers. 

The mi/ms ratio for the experiment can be determined by a linear combination of the results 
presented in figures and b. When sum rules are applied to the resultant composite map, it 
becomes possible to calculate the variation of mi/ms in reciprocal space, and two such maps 
are presented in figures and d. Surprisingly, despite the much lower strength of the EMCD 
signal originating from magnetite, these maps depend quite strongly on the assumed orbital angular 
momentum of its iron atoms. Figure]^ shows the variation of mi/ms ratio under the assumption of 
zero angular momentum on Fe atoms in the oxide. In this case, the mi/ms ratio varies very weakly 
within the area where measurements were performed, retaining the expected value of ~ 0.043. 
Stronger variations are observed only along the lines where the EMCD signal originating from bcc 
Fe is negligible, thus even tiny deviations lead to large changes. On the other hand, assuming that 
the iron atoms in the oxide layer have an enhanced orbital moment (for example mi = 1 /is), the 
picture changes substantially (figure [^). In the area close to the Thales circle (best represented 
by aperture position A) the value is significantly enhanced, while near the region best represented 
by aperture position B, the value remains close to the expected 0.043. The range of this variation 
scales with the size of orbital angular momentum on oxide iron atoms, with smaller values giving 
rise to less variation. Therefore, based on our simulations, the large variation of mi/ms observed in 
experiment suggests the presence of large, unquenched orbital angular momenta on the iron atoms 
in the oxide layer. 


Discussion 

The experimental design presented here provides two ways to study magnetic heterostructures that 
cannot be performed by any other method. First, by exploiting the angular dependency of the 
EMCD signal through the two different aperture pairings, it is possible to experimentally explore 
the different magnetic scattering contributions for both the metallic iron as well as its oxide surface 
layer. Second, by scanning the probe over a well-defined area from which individual datasets 
containing both spatial and spectral information are collected, it becomes possible to spatially 
segregate the EMCD signal with a spatial resolution of approximately 1.5 nm. 

By combining both of these methods, we reach the conclusion that the oxide layer is responsible 
for the measured changes of mi/ms- The simplest explanation could be that this is an artifact 
in the analysis due to the position and width of the interval used for post-edge normalization 
of the spectra. Particularly, based on discussion of the influence of magnetic EXAFS in XMCD 
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measurements of magnetite^, it was suggested that the normalization window should start above 
760 eV. We do observe that if we shrink the normalization interval to a much smaller width of lOeV, 
the calculated value of mi/ms becomes more sensitive to the position of normalization interval. 
However, note that the post-edge slope is very close to zero both aperture pairings over a wide 
energy range (see figure]^). As a result, as long as the normalization windows for both pairings 
are kept the same and their width is large enough to mitigate the influence of noise, pairing A will 
always have a larger value of mi/ms than pairing B. Hence a post-edge normalization artifact can 
be excluded as the reason for the enhanced mi/ms observed in aperture pairing A. 

Thus we are led explore the possibility that the orbital component to the net magnetization in 
the iron oxide layers is enhanced, as suggested by the simulations. Such results have some precedent 
in the literature. In the case of 7 -Fe 203 , Skoropata et ah have recently reported a strong increase 
in the orbital moment for iron in the outer shell of core-shell nanoparticles doped with cobalt^. For 
stoichiometric Fe 304 , Huang et ah reported orbital moments of 0.67 ± 0.07 hb at all temperatures 
measured using XMCD^i^, while Li et ah inferred a large orbital magnetic moment of 0.51 ±0.05 
at 10 K based on direct observation of the spin moment via Compton scattering and a comparison 
with literature magnetization data^. However, other studies suggest that the orbital moment is 
quenched within the stoichiometric Fe 304 , resulting in a nearly vanishing net valud®^, and these 
conflicting reports have lead to some controversy^^^®. 

Although our calculations provide strong support for the interpretation suggesting an enhanced 
orbital angular moment in oxide, we would like to present an alternative argument, which is inde¬ 
pendent of the simulations. Let’s assume that the orbital moments are actually compensated on 
both iron sublattices and in both iron valencies. The dynamical diffraction effect^^lP’ the 

contributions of the three sublattices in non-trivial way. A full disentangling of the individual con¬ 
tributions, as performed by Wang et. al.^^, is beyond the scope and intentions of this manuscript. 
However, in general, a linear combination of large spin moment contributions will yield a result 
larger than a corresponding linear combination of the small orbital moment contributions, espe¬ 
cially considering that there is a variation of the thickness of the oxide within the studied region. In 
such case, the observed mi/ms ratio would necessarily remain small. But this is in contradiction 
with our measurements, implying enhancement of the orbital moment in the oxide layer. On the 
other hand, it is important to point out that the exact value of mi/ms ratio for the oxide-dominated 
pairing A, or extracted from the oxide-rich area, is difficult to interpret without disentangling the 

individual site contribution^^lIlS! 

The reason for the variation of orbital magnetic moments in the literature is not always evident, 
although a number of theories exist. Kallmayer et ah used a monolayer sensitive XMCD technique 
to probe the interfaces of epitaxially-grown magnetite thin hlms on MgO and AI 2 O 3 substrates^. 
A modest enhancement of the orbital moment was observed at the interface between Fe 304 and 
AI 2 O 3 . They interpret this as a consequence of the reduced crystal symmetry of the magnetite at 
this interface arising from the incorporation of misfit dislocations due to the large lattice mismatch. 
Another explanation comes from an investigation into the potential for magnetite to harbor large, 
hidden orbital moments by E. Goering^. Goering calculates an average orbital moment of 1.0 
/xb per Fe atom within magnetite, but that it is nearly completely quenched. He concludes that 
“slight modihcations of the stoichiometry and crystallographic structure [of Magnetite] also give 
nonvanishing orbital moments.” Since the XMGD technique averages over very large areas of 
crystalline material, however, it is difficult to isolate these effects. 

Thus it appears that the observed enhancement of the orbital moment may be partially un¬ 
derstood by examining the literature. However, it is critical to emphasize that, due to the unique 
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nature of the EMCD technique presented here, it needs to be considered independently. For exam¬ 
ple, to the best of our knowledge, all of the quantitative magnetic information on stoichiometric 
Fe 304 to date comes from bulk systems or thin Elms where the signals are averaged over regions 
several hundreds of microns in diameter or more. The results presented here, on the other hand, are 
a summation of individual spectra acquired from volumes of material illuminated with an electron 
probe having a diameter of less than 1.5 nm. Thus nanoscale effects may play a greater role than 
for measurement techniques that probe much larger volumes of material. As an example of how 
this may manifest itself, we note that it is quite plausible that a series of correlated defects at the 
interface between Fe and magnetite could lead to a reduction in crystalline symmetry that may 
account for at least some enhancement of the orbital moment, and that the EMCD technique pre¬ 
sented here would be exceptionally sensitive to this. Moreover, although not explicitly investigated 
in this experiment, the preparation of the lamella with the FIB results in the presence of Gallium 
impurities on the exposed iron surface. This may influence the growth and magnetic behavior of the 
cladding oxide layers. For example. Gallium ions have been shown to be soluble in magnetite and 
take on a 3-|- oxidation state. They substitute for Fe^’*' at the tetrahedral sites, thereby leading to 
an enhancement of the net magnetic moment® Although the influence of these impurities on the 
orbital magnetic moment has not been directly investigated, the suppression of antiferrimagnetic 
balance of Fe^^ ions associated with their introduction may further contribute to the emergence of 
uncompensated orbital moments. 

In conclusion, we present an EMGD-based method that enables the quantitative analysis of 
magnetic moments in chemically and structurally distinct overlapping magnetic thin films in the 
TEM with nanoscale spatial resolution. Our structural and chemical analysis of the TEM lamella 
indicate that the system is best described as a free-standing trilayer of Fe 3_504 / Fe / Fe 3_504 with 
6 close to zero. By angularly and spatially segregating the EMGD signal, we observe significant 
variations in the measured value of mi/ms which can be correlated to the presence of these oxide 
layers. Simulations suggest that this can be understood if the orbital component of the net magnetic 
moment in the oxide is enhanced. Thus this method is capable of providing signihcant insight into 
the nature of nanoscale magnetism in a way not yet possible for any other technique. 


Methods 

Sample Preparation 

A thin him of pure iron was deposited by molecular beam epitaxy onto a single crystal (0 0 1) MgO 
substrate. The sample was prepared for the TEM using the FIB in-situ lift-out method®®. A 
protective bar of a Pt-G compound was locally deposited on the region of interest by decomposing 
a platinum carbon precursor gas in the presence of electrons accelerated to 3 kV. On top of this, a 
second bar of approximately 1 pm thickness was deposited by decomposing the same precursor gas 
in the presence of gallium ions accelerated to 30 kV. A lamella roughly 15 pm in length was milled, 
extracted, attached to a copper TEM grid using the precursor gas, and thinned to an estimated 
100 nm thickness using the 30 kV ion beam. A hnal polish using a 5 kV gallium ion beam was 
completed on both sides with an incidence angle of approximately 3° to the lamella face. 




Experimental Equipment 

The TEM and EELS measurements carried out as part of this study were conducted using a 
Tecnai F30 at Uppsala University as well as a Tecnai F20 at the Argonne National Laboratory 
Electron Microscopy Center (ANL EMCenter). Both instruments were equipped with Schottky 
Field Emission guns and operated at 300 kV and 200 kV, respectively. EELS data were acquired on 
a Tridiem Gatan Image Filter (GIF) (Gatan Inc.). The held emission gun was operated in such a 
way to produce a high current at the expense of energy resolution, which was close to 1.3 eV, taken 
as the Full Width Half Max (FWHM) value of the Zero Loss Peak (ZLP). All EELS spectra were 
acquired with the microscope set in Scanning TEM (STEM) mode with a calibrated camera length 
of 1877 mm. The minicondenser lens was switched off to yield a lower convergence angle of 1.6 mrad 
for a fully converged beam, which was diffraction limited by use of the secondary condenser aperture. 
We estimate that the spatial resolution of this conhguration is 1.2 nm. The EELS entrance aperture 
(physical diameter 1.0 mm) was used to set the collection angle at 2.4 mrad. The aperture was 
positioned using a script to excite the diffraction shift coils and the exact positions are shown in 
the insets of hgure |^,b. 

In STEM mode, a High Angle Annular Dark Field (HAADF) detector was used to produce a 
survey image using electrons scattered by angles larger than approximately 30 mrad. This resulted 
primarily in a mass-thickness contrast mechanism, but also included some diffraction contrast. The 
region from which all spectrum image data cubes were acquired is shown in green in hgure [T} A 
drift correction routine carried out at regular intervals ensured that the probe position within this 
region could be linked to the survey image and, subsequently, related between all of the individual 
data cubes. The probe was scanned perpendicular to the Fe/MgO interface with a pixel size of 
approximately 1.5 x 1.5 nm. 

The EMGD data were acquired by hrst shifting the aperture to the “Ghiral Plus” location of 
the pairing A shown in the inset of hgure and recording the energy range 480 - 890 eV using 
an acquisition time of 5 s per pixel and an energy dispersion of 0.2 eV per channel. The aperture 
position was subsequently shifted to the “Ghiral Minus” location and the same region was scanned 
again. This was repeated for both aperture positions in pairing B (see hgure S>. inset). For 
the core-loss EELS data, the dihraction pattern was shifted on-axis and the acquisition time was 
reduced to 2 s per pixel. The low-loss region was acquired with the dihraction pattern on-axis using 
a dispersion of 0.05 eV and an acquisition time of 1 x 10“^ s. Gritically, the survey image was 
not reacquired between aperture shifts, allowing for the same region to be scanned multiple times. 
Following the acquisition of each individual spectral image, correlated noise was accounted for by 
taking the average of 3\/N dark current measurements - where N is the total number of acquisitions 
- and subsequently subtracted from the gain normalized spectra^. 

Data Treatment 

Energy drift in the individual spectra within the spectral images was corrected for by using a 
cross-correlation algorithm to align the spectra to the Fe L 3 ionization edge within the regions of 
the him having the largest metallic iron content. The drift correction for the remaining spectra 
was determined by interpolating a spline ht between the non-corrected regions. For the EMGD 
spectra, the energy drift correction resulted in a gain averaging over approximately 15 channels, 
further improving the signal to noise ratio®'. Following this step, the ehects of plural scattering 
were removed by deconvolving all of the core-loss spectral images with the on-axis low-loss spectral 
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image^. The individual spectra were then summed into two datasets for each aperture position: 
the whole image and the “surface” region (see £gure[^. Following the summation, the pre-edge 
background was removed using a power-law background £t to the region between 670 - 700 e\®. 

The on-axis core-loss spectral image was used to quantify the relative amount of iron rpe and 
oxygen Tqx at any given pixel. This was accomplished by fitting the deconvolved data to the 
differential scattering cross-sections calculated using the Hartree-Slater model as implemented in 
Digital Micrograph. The effective electron mean free path was then computed for each pixel by 
utilizing the Fe:0 ratio R = rpe/rox at that position as proposed by Malis et al^ and implemented 
by Egerton in Matlab^S. This was used to calculate the absolute thickness tabs of the region at each 
pixel position by using the low-loss spectral image to hrst extract the relative thickness values. An 
equation governing the thickness of the oxide layer tox as a function of R and tabs was derived under 
the assumption of it being similar in structure to Fe 304 . The justihcation for this assumption is 
provided in the supplementary information. This yields equation 


Ax 


AbsPPe 


— D 
7 Pox 


7 Pox 


+ PFe 


( 1 ) 


where pox and ppe are the densities of the oxide and the metallic Fe layers, respectively, tpe was 
calculated as tabs — tox- The primary sources of systematic error for this calculation include the 
choice of cross-section model and the quantification routine, uncertainty in the local densities of 
the oxide and metal, the potential presence of an amorphous carbon coating layer, and unknown 
stoichiometric deviations from pure Fe 304 . Since many of these error sources are difficult to quantify, 
we are not able to provide systematic error bars for this calculation. Despite that, we note that the 
statistical error appears to be quite low due to the high signal to noise ratio. Consequently, we feel 
that these data provide a constructive qualitative assessment of the thicknesses of the individual 
metal and oxide layers. 

The EMCD signal was computed by hrst interpolating the background-removed data to a dis¬ 
persion of 0.01 eV and then aligning the two spectra of opposite chirality from each aperture pairing 
along the energy dispersion axis with a cross-correlation algorithm. The post-edge background was 
subsequently normalized to a window between 745 - 800 eV and the difference between the spectra 
was computed. With this, the conditions for sum rules are satished and rriL/ms can be calculated 
using equation 


rriL 

ms 


2 + J^,AI{E)dE 


( 2 ) 


3 AI{E)dE - 2 I(E)dE 

The mi/ms values in this paper were determined from the difference in the background-removed 
post-edge normalized data. The integration range for Fe L 3 was hxed at 700 - 715 eV. The lower 
bound of the integration range for Fe L 2 was set at 700 eV and the upper bound was varied from 
745 eV up to 800 eV. This variation is shown in figure |^. The reported values for mi/ms represent 
the average value of this variation between 745 - 800 eV while the error bars represent one standard 
deviation. 

We note that this way of processing data, particularly the post-edge normalization step before 
taking the difference, significantly suppresses the effects of asymmetry of the two beam case^^^^SP _ 
This follows from the cubic symmetry of both the magnetite and iron layer, which allows to write 
each ELNES spectrum as a linear combination of nonmagnetic spectrum N(E) and magnetic EMCD 
part Mz{E)\ 
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( 3 ) 


« A(n)N(E) + B(n)MJE) 

The post-edge normalization removes the differences in geometry-dependent coefficients ^4(^2), 
because EMCD is negligible in the post-edge region. Thus the difference of post-edge normalized 
spectra faithfully represents the EMCD spectrum. 

Simulations 

Simulations of the inelastic electron scattering were performed using the Bloch-waves method®^ 
utilizing the mats algorithmic. The orientations and thicknesses of the layers followed the observed 
ones. The transition matrix elements-mixed dynamical form factor^® were calculated using the 
operator maps techniqueC^ which allows for easy scaling of the magnitude of spin and orbital 
angular momentalC at the cost of being a dipole approximation. For bcc Fe, we used the values 
ms = l.QSps and = 0.086/1^ as reported by Chen et al.C For magnetite layers, we assumed 
spin moments of ms = 5.Ops and ms = 4.Ops, with a net magnetization of ms = 4.Ops- The 
value for mi was set to either zero (hgure [^) or Ips (hgure|^). The moments on the two iron 
sublattices of magnetite were oriented ferrimagnetically, assuming that its net magnetic moment is 
parallel with that of the iron layer. Resulting maps were scaled to take into account the thickness 
of each layer and amount of iron atoms per unit volume in iron (84.6 at./nm^) and magnetite (13.5 
and 26.99 at./nm^), respectively. Note the different magnitude of magnetic signals from the L 3 
edges in the iron (figure [^) and oxide (figure [^) layers, respectively. We have also performed test 
calculations comparing plane wave illumination and convergent beam illumination with a 1.6 mrad 
convergence angle beam in the (0 0 1) zone axis condition for bcc-Fe. This orientation was chosen 
so that the dynamical diffraction effects are maximized, thus providing an upper limit on the 
differences that could be expected between the two methods of calculation. However, thanks to 
Lorentzian broadening of the energy-hltered diffraction patterns due to energy loss proces^^, the 
resulting diffraction patterns and EMCD were very similar. Thus for the multilayer sample, we 
proceeded with simpler plane-wave simulations. 
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Figure 1: At top, an HAADF survey image denoting the region measured for EMCD is shown. The 
sub region denoted “surface” contains the largest volume percentage of oxide in the him and its 
position is marked as a red box on the maps. Relative composition maps for Fe and O extracted 
from this region are shown in the lower left. Absolute thickness maps for metallic Fe and Fe 3_504 
are shown in the lower right. Note that the values for the absolute thickness in the MgO substrate 
region are blacked out since they cannot be accurately calculated with this method. The pixel size 
of the maps is approximately 1.5 nm. 
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Figure 2: CBED pattern from the scanned region (green box in figure [^. Greyscale contrast is 
inverted to aid visualization in print. In (a), the raw pattern is shown. Low order indices for Fe 
(red) and Fe 304 (black) are provided in (b) and secondary reflections arising from the interference 
mechanism are colored yellow in (c). 
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Figure 3: (Color online) High resolution micrograph of the iron thin him in cross section. In (a), 
the HRTEM image is shown along with a box for the region from which an FFT is taken. The 
beam is oriented parallel to MgO [0 1 0] and Fe [1 10]. In (b), the FFT is shown as computed. In 
(c), the FFT is indexed using the hrst order spatial frequencies for Fe and Fe 304 . The unindexed 
spatial frequencies can be accounted for through convolution of the two materials, as discussed in 
the supplementary information. This information allows us to construct a structural model for the 
sample in the TEM, which is presented in (d). 
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Figure 4: Background subtracted and post-edge normalized EELS spectra for both chiral locations 
in aperture pairings A (a) and B (b). Both spectra are normalized to the max value of the Chiral 
Plus spectra and the normalization factor is provided. The difference spectrum (EMCD signal) is 
also presented. Inset is a schematic depicting the positions of the different aperture pairings in 
the reciprocal plane. The individual aperture positions are marked -|- and - to denote Chiral Plus 
and Chiral Minus, respectively, (c) Integral of the EMCD signal for both aperture positions, (d) 
Variation of mi/ms as a function of integration range. 
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Figure 5: Background subtracted and post-edge normalized EELS spectra for the summation over 
the “surface” region, as denoted in hgure 1. Aperture pairings A and B are presented in (a) and 
(b), respectively. The EMCD signals for each chiral pair are also presented. 
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Figure 6 : Top row: calculated maps of the EMCD signal at the L 3 edge of iron originating from 
bcc Fe (a) and two 12.5 nm thick oxide layers (b). Bottom row: maps showing the variation of the 
resulting mi/ms ratio, assuming zero orbital angular momentum in oxide (c) or mi = l/i^ (d). 
The Thales circle is overlaid on all maps. 
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